Machine learning approaches delimit cryptic taxa in a previously intractable species complex.

Cryptic species are not diagnosable via morphological criteria, but can be detected through analysis of DNA sequences. A number of methods have been developed for identifying species based on genetic data; however, these methods are prone to over-splitting taxa with extreme population structure, such as dispersal-limited organisms. Machine learning methodologies have the potential to overcome this challenge. Here, we apply such approaches, using a large dataset generated through hybrid target enrichment of ultraconserved elements (UCEs). Our study taxon is the Aoraki denticulata species complex, a lineage of extremely low-dispersal arachnids endemic to the South Island of Aotearoa New Zealand. This group of mite harvesters has been the subject of previous species delimitation studies using smaller datasets generated through Sanger sequencing and analytical approaches that rely on multispecies coalescent models and barcoding gap discovery. Those analyses yielded a number of putative cryptic species that seems unrealistic and extreme, based on what we know about species' geographic ranges and genetic diversity in non-cryptic mite harvesters. We find that machine learning approaches, on the other hand, identify cryptic species with geographic ranges that are similar to those seen in other morphologically diagnosable mite harvesters in Aotearoa New Zealand's South Island. We performed both unsupervised and supervised machine learning analyses, the latter with training data drawn either from animals broadly (vagile and non-vagile) or from a custom training dataset from dispersal-limited harvesters. We conclude that applying machine learning approaches to the analysis of UCE-derived genetic data is an effective method for delimiting species in complexes of low-vagility cryptic species, and that the incorporation of training data from biologically relevant analogues can be critically informative.

Under the hood: Phylogenomics of hooded tick spiders (Arachnida, Ricinulei) uncovers discordance between morphology and molecules.

Ricinulei or hooded tick-spiders are a cryptic and ancient group of arachnids. The order consists of around 100 highly endemic extant species restricted to the Afrotropics and the Neotropics along with 22 fossil species. Their antiquity and low vagility make them an excellent group with which to interrogate biogeographic questions. To date, only four molecular analyses have been conducted on the group and they failed to resolve the relationships of the main lineages and even recovering the non-monophyly of the three genera. These studies were limited to a few Sanger loci or phylogenomic analyses with at most seven ingroup samples. To increase phylogenetic resolution in this little-understood and poorly studied group, we present the most comprehensive phylogenomic study of Ricinulei to date leveraging the Arachnida ultra-conserved element probe set. With a data set of 473 loci across 96 ingroup samples, analyses resolved a monophyletic Neotropical clade consisting of four main lineages. Two of them correspond to the current genera Cryptocellus and Pseudocellus while topology testing revealed one lineage to likely be a phylogenetic reconstruction artefact. The fourth lineage, restricted to Northwestern, Andean South America, is consistent with the Cryptocellus magnus group, likely corresponding to the historical genus Heteroricinoides. Since we did not sample the type species for this old genus, we do not formally re-erect Heteroricinoides but our data suggest the need for a thorough morphological re-examination of Neotropical Ricinulei.

Dynamic evolution of locomotor performance independent of changes in extended phenotype use in spiders.

Many animals use self-built structures (extended phenotypes) to enhance body functions, such as thermoregulation, prey capture or defence. Yet, it is unclear whether the evolution of animal constructions supplements or substitutes body functions-with disparate feedbacks on trait evolution. Here, using brown spiders (Araneae: marronoid clade), we explored if the evolutionary loss and gain of silken webs as extended prey capture devices correlates with alterations in traits known to play an important role in predatory strikes-locomotor performance (sprint speed) and leg spination (expression of capture spines on front legs). We found that in this group high locomotor performance, with running speeds of over 100 body lengths per second, evolved repeatedly-both in web-building and cursorial spiders. There was no correlation with running speed, and leg spination only poorly correlated, relative to the use of extended phenotypes, indicating that web use does not reduce selective pressures on body functions involved in prey capture and defence per se. Consequently, extended prey capture devices serve as supplements rather than substitutions to body traits and may only be beneficial in conjunction with certain life-history traits, possibly explaining the rare evolution and repeated loss of trapping strategies in predatory animals.

An Opiliones-specific ultraconserved element probe set with a near-complete family-level phylogeny.

Sequence capture of ultraconserved elements (UCEs) has transformed molecular systematics across many taxa, with arachnids being no exception. The probe set available for Arachnida has been repeatedly used across multiple arachnid lineages and taxonomic levels, however more specific probe sets for spiders have demonstrated that more UCEs can be recovered with higher probe specificity. In this study, we develop an Opiliones-specific UCE probe set targeting 1915 UCEs using a combination of probes designed from genomes and transcriptomes, as well as the most useful probes from the Arachnida probe set. We demonstrate the effectiveness of this probe set across Opiliones with the most complete family-level phylogeny made to date, including representatives from 61 of 63 currently described families. We also test UCE recovery from historical specimens with degraded DNA, examine population-level data sets, and assess "backwards compatibility" with samples hybridized with the Arachnida probe set. The resulting phylogenies - which include specimens hybridized using both the Opiliones and Arachnida probe sets, historical specimens, and transcriptomes - are largely congruent with previous multi-locus and phylogenomic analyses. The probe set is also "backwards compatible", increasing the number of loci obtained in samples previously hybridized with the Arachnida probe set, and shows high utility down to shallow population-level divergences. This probe set has the potential to further transform Opiliones molecular systematics, resolving many long-standing taxonomic issues plaguing this lineage.

Using natural history to guide supervised machine learning for cryptic species delimitation with genetic data.

The diversity of biological and ecological characteristics of organisms, and the underlying genetic patterns and processes of speciation, makes the development of universally applicable genetic species delimitation methods challenging. Many approaches, like those incorporating the multispecies coalescent, sometimes delimit populations and overestimate species numbers. This issue is exacerbated in taxa with inherently high population structure due to low dispersal ability, and in cryptic species resulting from nonecological speciation. These taxa present a conundrum when delimiting species: analyses rely heavily, if not entirely, on genetic data which over split species, while other lines of evidence lump. We showcase this conundrum in the harvester Theromaster brunneus, a low dispersal taxon with a wide geographic distribution and high potential for cryptic species. Integrating morphology, mitochondrial, and sub-genomic (double-digest RADSeq and ultraconserved elements) data, we find high discordance across analyses and data types in the number of inferred species, with further evidence that multispecies coalescent approaches over split. We demonstrate the power of a supervised machine learning approach in effectively delimiting cryptic species by creating a "custom" training data set derived from a well-studied lineage with similar biological characteristics as Theromaster. This novel approach uses known taxa with particular biological characteristics to inform unknown taxa with similar characteristics, using modern computational tools ideally suited for species delimitation. The approach also considers the natural history of organisms to make more biologically informed species delimitation decisions, and in principle is broadly applicable for taxa across the tree of life.

First record of Travunioidea (Arachnida: Opiliones: Laniatores) from China, with the description of a new monotypic genus from a cave.

A new genus and species of troglobitic Opiliones, Sinonychia martensi gen. n. sp. n., is described and illustrated. This is the first report of the superfamily Travunioidea Absolon Kratochvíl, 1932 from China. The new species is assigned to Cladonychiidae Hadzi, 1935 on the basis of morphology of midgut and male genitalia, and can be easily distinguished from related genera by several characters.

An approach using ddRADseq and machine learning for understanding speciation in Antarctic Antarctophilinidae gastropods.

Sampling impediments and paucity of suitable material for molecular analyses have precluded the study of speciation and radiation of deep-sea species in Antarctica. We analyzed barcodes together with genome-wide single nucleotide polymorphisms obtained from double digestion restriction site-associated DNA sequencing (ddRADseq) for species in the family Antarctophilinidae. We also reevaluated the fossil record associated with this taxon to provide further insights into the origin of the group. Novel approaches to identify distinctive genetic lineages, including unsupervised machine learning variational autoencoder plots, were used to establish species hypothesis frameworks. In this sense, three undescribed species and a complex of cryptic species were identified, suggesting allopatric speciation connected to geographic or bathymetric isolation. We further observed that the shallow waters around the Scotia Arc and on the continental shelf in the Weddell Sea present high endemism and diversity. In contrast, likely due to the glacial pressure during the Cenozoic, a deep-sea group with fewer species emerged expanding over great areas in the South-Atlantic Antarctic Ridge. Our study agrees on how diachronic paleoclimatic and current environmental factors shaped Antarctic communities both at the shallow and deep-sea levels, promoting Antarctica as the center of origin for numerous taxa such as gastropod mollusks.

Sequence Capture Phylogenomics of True Spiders Reveals Convergent Evolution of Respiratory Systems.

The common ancestor of spiders likely used silk to line burrows or make simple webs, with specialized spinning organs and aerial webs originating with the evolution of the megadiverse "true spiders" (Araneomorphae). The base of the araneomorph tree also concentrates the greatest number of changes in respiratory structures, a character system whose evolution is still poorly understood, and that might be related to the evolution of silk glands. Emphasizing a dense sampling of multiple araneomorph lineages where tracheal systems likely originated, we gathered genomic-scale data and reconstructed a phylogeny of true spiders. This robust phylogenomic framework was used to conduct maximum likelihood and Bayesian character evolution analyses for respiratory systems, silk glands, and aerial webs, based on a combination of original and published data. Our results indicate that in true spiders, posterior book lungs were transformed into morphologically similar tracheal systems six times independently, after the evolution of novel silk gland systems and the origin of aerial webs. From these comparative data, we put forth a novel hypothesis that early-diverging web-building spiders were faced with new energetic demands for spinning, which prompted the evolution of similar tracheal systems via convergence; we also propose tests of predictions derived from this hypothesis.[Book lungs; discrete character evolution; respiratory systems; silk; spider web evolution; ultraconserved elements.].

Pleistocene persistence and expansion in tarantulas on the Colorado Plateau and the effects of missing data on phylogeographical inferences from RADseq.

Montane species endemic to the "sky islands" of the North American southwest were significantly impacted by changing climates during the Pleistocene. We combined mitochondrial and genomic data with species distribution modelling to determine whether Aphonopelma marxi, a large tarantula from the nearby Colorado Plateau, was similarly impacted by glacial climates. Genetic analyses revealed that the species comprises three main clades that diverged in the Pleistocene. A clade distributed along the Mogollon Rim appears to have persisted in place during glacial conditions, whereas the other two clades probably colonized central and northeastern portions of the species' range from refugia in canyons. Climate models support this hypothesis for the Mogollon Rim, but late glacial climate data appear too coarse to detect suitable areas in canyons. Locations of canyon refugia could not be inferred from genomic analyses due to missing data, encouraging us to explore the effect of missing loci in phylogeographical inferences using RADseq. Results from analyses with varying amounts of missing data suggest that samples with large amounts of missing data can still improve inferences, and the specific loci that are missing matters more than the number of missing loci. This study highlights the profound impact of Pleistocene climates on tarantulas endemic to the Colorado Plateau, as well as the mixed nature of the region's fauna. Some animals recently colonized from nearby deserts as glacial climates receded, whereas others, like tarantulas, appear to have persisted on the Mogollon Rim and in refugia associated with the region's famous river-cut canyons.

Integrative species delimitation reveals cryptic diversity in the southern Appalachian Antrodiaetus unicolor (Araneae: Antrodiaetidae) species complex.

Although species delimitation can be highly contentious, the development of reliable methods to accurately ascertain species boundaries is an imperative step in cataloguing and describing Earth's quickly disappearing biodiversity. Spider species delimitation remains largely based on morphological characters; however, many mygalomorph spider populations are morphologically indistinguishable from each other yet have considerable molecular divergence. The focus of our study, the Antrodiaetus unicolor species complex containing two sympatric species, exhibits this pattern of relative morphological stasis with considerable genetic divergence across its distribution. A past study using two molecular markers, COI and 28S, revealed that A. unicolor is paraphyletic with respect to A. microunicolor. To better investigate species boundaries in the complex, we implement the cohesion species concept and use multiple lines of evidence for testing genetic exchangeability and ecological interchangeability. Our integrative approach includes extensively sampling homologous loci across the genome using a RADseq approach (3RAD), assessing population structure across their geographic range using multiple genetic clustering analyses that include structure, principal components analysis and a recently developed unsupervised machine learning approach (Variational Autoencoder). We evaluate ecological similarity by using large-scale ecological data for niche-based distribution modelling. Based on our analyses, we conclude that this complex has at least one additional species as well as confirm species delimitations based on previous less comprehensive approaches. Our study demonstrates the efficacy of genomic-scale data for recognizing cryptic species, suggesting that species delimitation with one data type, whether one mitochondrial gene or morphology, may underestimate true species diversity in morphologically homogenous taxa with low vagility.

Sitticine jumping spiders: phylogeny, classification, and chromosomes (Araneae, Salticidae, Sitticini).

The systematics of sitticine jumping spiders is reviewed, with a focus on the Palearctic and Nearctic regions, in order to revise their generic classification, clarify the species of one region (Canada), and study their chromosomes. A genome-wide molecular phylogeny of 23 sitticine species, using more than 700 loci from the arachnid Ultra-Conserved Element (UCE) probeset, confirms the Neotropical origins of sitticines, whose basal divergence separates the new subtribeAillutticina (a group of five Neotropical genera) from the subtribe Sitticina (five genera of Eurasia and the Americas). The phylogeny shows that most Eurasian sitticines form a relatively recent and rapid radiation, which we unite into the genus Attulus Simon, 1868, consisting of the subgenera Sitticus Simon, 1901 (seven described species), Attulus (41 described species), and Sittilong Prószynski, 2017 (one species). Five species of Attulus occur natively in North America, presumably through dispersals back from the Eurasian radiation, but an additional three species were more recently introduced from Eurasia. Attus palustris Peckham & Peckham, 1883 is considered to be a full synonym of Euophrys floricola C. L. Koch, 1837 (not a distinct subspecies). Attus sylvestris Emerton, 1891 is removed from synonymy and recognized as a senior synonym of Sitticus magnus Chamberlin & Ivie, 1944. Thus, the five native Attulus in North America are Attulus floricola, A. sylvestris, A. cutleri, A. striatus, and A. finschi. The other sitticines of Canada and the U.S.A. are placed in separate genera, all of which arose from a Neotropical radiation including Jollas Simon, 1901 and Tomis F.O.Pickard-Cambridge, 1901: (1) Attinella Banks, 1905 (A. dorsata, A. concolor, A. juniperi), (2) Tomis (T. welchi), and (3) Sittisax Prószynski, 2017 (S. ranieri). All Neotropical and Caribbean "Sitticus" are transferred to either Jollas (12 species total) or Tomis (14 species). Attinella (three species) and Tomis are both removed from synonymy with Sitticus; the synonymy of Sitticus cabellensis Prószynski, 1971 with Pseudattulus kratochvili Caporiacco, 1947 is restored; Pseudattulus Caporiacco, 1947 is synonymized with Tomis. Six generic names are newly synonymized with Attulus and one with Attinella. Two Neotropical species are described as new, Jollas cupreus sp. nov. and Tomis manabita sp. nov. Forty-six new combinations are established and three are restored. Three species synonymies are restored, one is new, and two are rejected. Across this diversity of species is a striking diversification of chromosome complements, with X-autosome fusions occurring at least four times to produce neo-Y sex chromosome systems (X1X2Y and X1X2X3Y), some of which (Sittisax ranieri and S. saxicola) are sufficiently derived as to no longer preserve the simple traces of ancestral X material. The correlated distribution of neo-Y and a base autosome number of 28 suggests that neo-Y origins occurred preferentially in lineages with the presence of an extra pair of autosomes.

Sequence capture phylogenomics of historical ethanol-preserved museum specimens: Unlocking the rest of the vault.

Natural history collections play a crucial role in biodiversity research, and museum specimens are increasingly being incorporated into modern genetics-based studies. Sequence capture methods have proven incredibly useful for phylogenomics, providing the additional ability to sequence historical museum specimens with highly degraded DNA, which until recently have been deemed less valuable for genetic work. The successful sequencing of ultraconserved elements (UCEs) from historical museum specimens has been demonstrated on multiple tissue types including dried bird skins, formalin-fixed squamates and pinned insects. However, no study has thoroughly demonstrated this approach for historical ethanol-preserved museum specimens. Alongside sequencing of "fresh" specimens preserved in >95% ethanol and stored at -80°C, we used extraction techniques specifically designed for degraded DNA coupled with sequence capture protocols to sequence UCEs from historical museum specimens preserved in 70%-80% ethanol and stored at room temperature, the standard for such ethanol-preserved museum collections. Across 35 fresh and 15 historical museum samples of the arachnid order Opiliones, an average of 345 UCE loci were included in phylogenomic matrices, with museum samples ranging from six to 495 loci. We successfully demonstrate the inclusion of historical ethanol-preserved museum specimens in modern sequence capture phylogenomic studies, show a high frequency of variant bases at the species and population levels, and from off-target reads successfully recover multiple loci traditionally sequenced in multilocus studies including mitochondrial loci and nuclear rRNA loci. The methods detailed in this study will allow researchers to potentially acquire genetic data from millions of ethanol-preserved museum specimens held in collections worldwide.

A demonstration of unsupervised machine learning in species delimitation.

One major challenge to delimiting species with genetic data is successfully differentiating population structure from species-level divergence, an issue exacerbated in taxa inhabiting naturally fragmented habitats. Many fields of science are now using machine learning, and in evolutionary biology supervised machine learning has recently been used to infer species boundaries. These supervised methods require training data with associated labels. Conversely, unsupervised machine learning (UML) uses inherent data structure and does not require user-specified training labels, potentially providing more objectivity in species delimitation. In the context of integrative taxonomy, we demonstrate the utility of three UML approaches (random forests, variational autoencoders, t-distributed stochastic neighbor embedding) for species delimitation in an arachnid taxon with high population genetic structure (Opiliones, Laniatores, Metanonychus). We find that UML approaches successfully cluster samples according to species-level divergences and not high levels of population structure, while model-based validation methods severely over-split putative species. UML offers intuitive data visualization in two-dimensional space, the ability to accommodate various data types, and has potential in many areas of systematic and evolutionary biology. We argue that machine learning methods are ideally suited for species delimitation and may perform well in many natural systems and across taxa with diverse biological characteristics.

Phylogenomic analysis and revised classification of atypoid mygalomorph spiders (Araneae, Mygalomorphae), with notes on arachnid ultraconserved element loci.

The atypoid mygalomorphs include spiders from three described families that build a diverse array of entrance web constructs, including funnel-and-sheet webs, purse webs, trapdoors, turrets and silken collars. Molecular phylogenetic analyses have generally supported the monophyly of Atypoidea, but prior studies have not sampled all relevant taxa. Here we generated a dataset of ultraconserved element loci for all described atypoid genera, including taxa (Mecicobothrium and Hexurella) key to understanding familial monophyly, divergence times, and patterns of entrance web evolution. We show that the conserved regions of the arachnid UCE probe set target exons, such that it should be possible to combine UCE and transcriptome datasets in arachnids. We also show that different UCE probes sometimes target the same protein, and under the matching parameters used here show that UCE alignments sometimes include non-orthologs. Using multiple curated phylogenomic matrices we recover a monophyletic Atypoidea, and reveal that the family Mecicobothriidae comprises four separate and divergent lineages. Fossil-calibrated divergence time analyses suggest ancient Triassic (or older) origins for several relictual atypoid lineages, with late Cretaceous/early Tertiary divergences within some genera indicating a high potential for cryptic species diversity. The ancestral entrance web construct for atypoids, and all mygalomorphs, is reconstructed as a funnel-and-sheet web.

Sequence capture phylogenomics of eyeless Cicurina spiders from Texas caves, with emphasis on US federally-endangered species from Bexar County (Araneae, Hahniidae).

Morphological, mitochondrial, and nuclear phylogenomic data were combined to address phylogenetic and species delimitation questions in cave-limited Cicurina spiders from central Texas. Special effort was focused on specimens and cave locations in the San Antonio region (Bexar County), home to four eyeless species listed as US Federally Endangered. Sequence capture experiments resulted in the recovery of ~200-400 homologous ultra-conserved element (UCE) nuclear loci across taxa, and nearly complete COI mitochondrial DNA sequences from the same set of individuals. Some of these nuclear and mitochondrial sequences were recovered from "standard" museum specimens without special preservation of DNA material, including museum specimens preserved in the 1990s. Multiple phylogenetic analyses of the UCE data agree in the recovery of two major lineages of eyeless Cicurina in Texas. These lineages also differ in mitochondrial clade membership, female genitalic morphology, degree of troglomorphy (as measured by relative leg length), and are mostly allopatric across much of Texas. Rare sympatry was confirmed in Bexar County, where members of the two major clades sometimes co-exist in the same karst feature. Both nuclear phylogenomic and mitochondrial data indicate the existence of undescribed species from the San Antonio region, although further sampling and collection of adult specimens is needed to explicitly test these hypotheses. Our data support the two following species synonymies (Cicurina venii Gertsch, 1992 = Cicurina madla Gertsch, 1992; Cicurina loftini Cokendolpher, 2004 = Cicurina vespera Gertsch, 1992), formally proposed here. Overall, our taxonomy-focused research has many important conservation implications, and again highlights the fundamental importance of robust taxonomy in conservation research.

A stable phylogenomic classification of Travunioidea (Arachnida, Opiliones, Laniatores) based on sequence capture of ultraconserved elements.

Molecular phylogenetics has transitioned into the phylogenomic era, with data derived from next-generation sequencing technologies allowing unprecedented phylogenetic resolution in all animal groups, including understudied invertebrate taxa. Within the most diverse harvestmen suborder, Laniatores, most relationships at all taxonomic levels have yet to be explored from a phylogenomics perspective. Travunioidea is an early-diverging lineage of laniatorean harvestmen with a Laurasian distribution, with species distributed in eastern Asia, eastern and western North America, and south-central Europe. This clade has had a challenging taxonomic history, but the current classification consists of ~77 species in three families, the Travuniidae, Paranonychidae, and Nippononychidae. Travunioidea classification has traditionally been based on structure of the tarsal claws of the hind legs. However, it is now clear that tarsal claw structure is a poor taxonomic character due to homoplasy at all taxonomic levels. Here, we utilize DNA sequences derived from capture of ultraconserved elements (UCEs) to reconstruct travunioid relationships. Data matrices consisting of 317-677 loci were used in maximum likelihood, Bayesian, and species tree analyses. Resulting phylogenies recover four consistent and highly supported clades; the phylogenetic position and taxonomic status of the enigmatic genus Yuria is less certain. Based on the resulting phylogenies, a revision of Travunioidea is proposed, now consisting of the Travuniidae, Cladonychiidae, Paranonychidae (Nippononychidae is synonymized), and the new family Cryptomastridae Derkarabetian & Hedin, fam. n., diagnosed here. The phylogenetic utility and diagnostic features of the intestinal complex and male genitalia are discussed in light of phylogenomic results, and the inappropriateness of the tarsal claw in diagnosing higher-level taxa is further corroborated.

Phylogenomic reclassification of the world's most venomous spiders (Mygalomorphae, Atracinae), with implications for venom evolution.

Here we show that the most venomous spiders in the world are phylogenetically misplaced. Australian atracine spiders (family Hexathelidae), including the notorious Sydney funnel-web spider Atrax robustus, produce venom peptides that can kill people. Intriguingly, eastern Australian mouse spiders (family Actinopodidae) are also medically dangerous, possessing venom peptides strikingly similar to Atrax hexatoxins. Based on the standing morphology-based classification, mouse spiders are hypothesized distant relatives of atracines, having diverged over 200 million years ago. Using sequence-capture phylogenomics, we instead show convincingly that hexathelids are non-monophyletic, and that atracines are sister to actinopodids. Three new mygalomorph lineages are elevated to the family level, and a revised circumscription of Hexathelidae is presented. Re-writing this phylogenetic story has major implications for how we study venom evolution in these spiders, and potentially genuine consequences for antivenom development and bite treatment research. More generally, our research provides a textbook example of the applied importance of modern phylogenomic research.

Cryptic elevational zonation in trapdoor spiders (Araneae, Antrodiaetidae, Aliatypus janus complex) from the California southern Sierra Nevada.

The relative roles of ecological niche conservatism versus niche divergence in promoting montane speciation remains an important topic in biogeography. Here, our aim was to test whether lineage diversification in a species complex of trapdoor spiders corresponds with riverine barriers or with an ecological gradient associated with elevational tiering. Aliatypus janus was sampled from throughout its range, with emphasis on populations in the southern Sierra Nevada Mountains of California. We collected multi-locus genetic data to generate a species tree for A. janus and its close relatives. Coalescent based hypothesis tests were conducted to determine if genetic breaks within A. janus conform to riverine barriers. Ecological niche models (ENM) under current and Last Glacial Maximum (LGM) conditions were generated and hypothesis tests of niche conservatism and divergence were performed. Coalescent analyses reveal deeply divergent genetic lineages within A. janus, likely corresponding to cryptic species. Two primary lineages meet along an elevational gradient on the western slopes of the southern Sierra Nevada Mountains. ENMs under both current and LGM conditions indicate that these groups occupy largely non-overlapping niches. ENM hypothesis testing rejected niche identity between the two groups, and supported a sharp ecological gradient occurring where the groups meet. However, the niche similarity test indicated that the two groups may not inhabit different background niches. The Sierra Nevada Mountains provide a natural laboratory for simultaneously testing ecological niche divergence and conservatism and their role in speciation across a diverse range of taxa. Aliatypus janus represents a species complex with cryptic lineages that may have diverged due to parapatric speciation along an ecological gradient, or been maintained by the evolution of ecological niche differences following allopatric speciation.

Comparative performance of double-digest RAD sequencing across divergent arachnid lineages.

Next-generation sequencing technologies now allow researchers of non-model systems to perform genome-based studies without the requirement of a (often unavailable) closely related genomic reference. We evaluated the role of restriction endonuclease (RE) selection in double-digest restriction-site-associated DNA sequencing (ddRADseq) by generating reduced representation genome-wide data using four different RE combinations. Our expectation was that RE selections targeting longer, more complex restriction sites would recover fewer loci than RE with shorter, less complex sites. We sequenced a diverse sample of non-model arachnids, including five congeneric pairs of harvestmen (Opiliones) and four pairs of spiders (Araneae). Sample pairs consisted of either conspecifics or closely related congeneric taxa, and in total 26 sample pair analyses were tested. Sequence demultiplexing, read clustering and variant calling were performed in the pyRAD program. The 6-base pair cutter EcoRI combined with methylated site-specific 4-base pair cutter MspI produced, on average, the greatest numbers of intra-individual loci and shared loci per sample pair. As expected, the number of shared loci recovered for a sample pair covaried with the degree of genetic divergence, estimated with cytochrome oxidase I sequences, although this relationship was non-linear. Our comparative results will prove useful in guiding protocol selection for ddRADseq experiments on many arachnid taxa where reference genomes, even from closely related species, are unavailable.

High phylogenetic utility of an ultraconserved element probe set designed for Arachnida.

Arachnida is an ancient, diverse and ecologically important animal group that contains a number of species of interest for medical, agricultural and engineering applications. Despite their importance, many aspects of the arachnid tree of life remain unresolved, hindering comparative approaches to arachnid biology. Biologists have made considerable efforts to resolve the arachnid phylogeny; yet, limited and challenging morphological characters, as well as a dearth of genetic resources, have hindered progress. Here, we present a genomic toolkit for arachnids featuring hundreds of conserved DNA regions (ultraconserved elements or UCEs) that allow targeted sequencing of any species in the arachnid tree of life. We used recently developed capture probes designed from conserved regions of available arachnid genomes to enrich a sample of loci from 32 diverse arachnids. Sequence capture returned an average of 487 UCE loci for all species, with a range from 170 to 722. Phylogenetic analysis of these UCEs produced a highly resolved arachnid tree with relationships largely consistent with recent transcriptome-based phylogenies. We also tested the phylogenetic informativeness of UCE probes within the spider, scorpion and harvestman orders, demonstrating the utility of these markers at shallower taxonomic scales and suggesting that these loci will be useful for species-level differences. This probe set will open the door to phylogenomic and population genomic studies across the arachnid tree of life, enabling systematics, species delimitation, species discovery and conservation of these diverse arthropods.